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Abstract 

The orbital magnetization of tlie electron gas on a two-dimensional kagome lattice under a per- 
pendicular magnetic field is theoretically investigated. The interplay between the lattice geometry 
and magnetic field induce nontrivial A;-space Chern invariant in the magnetic Brillouin zone, which 
turns to result in profound effects on the magnetization properties. We show that the Berry-phase 
term in the magnetization gives a paramagnetic contribution, while the conventional term brought 
about by the magnetic response of the magnetic Bloch bands produces a diamagnetic contribution. 
As a result, the superposition of these two components gives rise to a delicate oscillatory struc- 
ture in the magnetization curve when varying the electron filling factor. The relationship between 
this oscillatory behavior and the Hofstadter energy spectrum is revealed by selectively discussing 
the magnetization and its two components at the commensurate fluxes of /=l/4, 1/3, and 1/6, 
respectively. In particular, we reveal as a typical example the fractal structure in the magnetic 
oscillations by tuning the commensurate flux around /=l/4. The finite-temperature effect on the 
magnetization is also discussed. 
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I. INTRODUCTION 



Over the past three decades, the orbital dynamics of two-dimensional (2D) electrons 
coupled to a uniform perpendicular magnetic field has been a subject of special interest 
in condensed matter physics. Great attention has been paid to investigate the well-known 
Hofstadter butterfly-like diagram, a remarkable intricate, detailed, and self-similar feature 
for the subband spectrum of single-particle eigenenergies, in different 2D structures such 
as square triangular honeycomb jsl, and kagome lattices {4], lateral superlattices 
{st, superstructures with fiat bands jst, quasiperiodic tilings f^], fractal networks js], and 
rhombus tilings jol- This magnetic- field-induced frustration in 2D structures, caused by the 
competition between the periodic potential imposed by the lattice structure and the length 
scale provided by the magnetic field, is the headstream of rich and novel physical phenomena 
in these 2D systems. Among them the orbital magnetization is a very interesting one and 
needs to be paid special attention due to its intrinsic relationship with the unique topological 
structure brought about by the interplay between the special lattice geometry and magnetic 
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subbands. Unfortunately, to date, there are very scarce studies 
devote to enlightening this revealing relationship. 

Motivated by this observation, in this paper, we study the orbital magnetization proper- 
ties of the electron gas on a two-dimensional kagome lattice under a perpendicular magnetic 
field. Since our attention is solely on the orbital character in magnetization brought about 
by the interplay between the kagome lattice structure and the magnetic field, thus unlike 
most of previous work, the kagome lattice considered in this paper is spin nonmagnetic in 
itself, namely, no spin structure is exposed in the absence of the external magnetic field. 
The nonmagnetic kagome lattice structure has been either fabricated by modern patterning 
techniques II2I . fisj ] or observed in reconstructed semiconductor surfaces [14]] . In the former 
case, remarkably, the electron filling factor (namely, the Fermi energy) can be readily con- 
trolled by applying a gate voltage ISj. Our lattice model is free from the constraint imposed 
on the k-p approximation used in the extensively studied GaAs two-dimensional electron 
gas (2DEG), in which the k-p Hamiltonian is only valid around the F point in the Brillouin 
zone (BZ). In contrast, our lattice model allows for any electron filling, which result in var- 
ious Fermi-surface topologies in the magnetic BZ, which in turn, as will be shown below, 
produces profound effects on the orbital magnetization properties and the related transport 
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The orbital magnetization for Bloch electrons combines two terms 
conventional orbital magnetic moment to characterize the magnetic response of the single- 
particle energy spectrum, and the other is a Berry-phase correction. These two terms play 
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. By varying the Fermi energy, we 



different roles in metallic and insulating regions 
obtain the follows: (i) Both the magnitudes of these two terms increase (decrease) with 
increasing (decreasing) the Fermi energy in metallic regions; (ii) In insulating regions the 
conventional term keeps a constant unchanged value (quantum step), while the Berry- phase 
linearly varies with the Fermi energy with a slope proportional to the system's Hall conduc- 
tance; (iii) A general fractal structure may occur in the orbital magnetization curve (as a 
function of the Fermi energy) by tuning the commensurate flux. 

The rest of this paper is organized as follows. In Sec. II we give the Hofstadter spectrum 
of the kagome lattice by diagonalizing the corresponding tight-binding Hamiltonian in an 
external magnetic field. In Sec. Ill the orbital magnetization and its two components are 
systematically discussed by varying various system parameters. A summary is given in Sec. 
IV. 



II. KAGOME LATTICE IN A MAGNETIC FIELD 



The kagome lattice is a two-dimensional periodic array of corner- sharing triangles with 
three sites per unit cell, as illustrated in Fig. [H Here a is the triangle edge length (we set 
it as the length unit in the following), and A, B, and C denote the three sites in a unit cell. 
Let us begin with the tight-binding Hamiltonian for the spinless electrons on the kagome 
lattice with zero magnetic field (B=0), which is given by 

n = tJ2clc,, (1) 

where t is hopping amplitude between the nearest-neighbor link and cj (q) is the cre- 
ation (annihilation) operator of an electron on lattice i. Hamiltonian ([1]) can be diagonalized 
in the momentum space as 

H = t^vl/tiJ(k)vl/k, (2) 

k 
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where \l/k=(ckA, c^b, c^c)^ is the three-component electron field operator. Each component 
of is the Fourier transform of q, i.e., 



mn 



where we have changed notation i — )■ (mns) by using (mn) to label the kagome unit cells 
and s to label the three sites in a unit cell. H{k) is given by 



/ 1 A 

' pI pi^ 



(4) 



pI pI 
\pi pI J 

where p]^=2 cos (k ■ sli) and ai=(— 1/2, — v^/2), a2=(l,0), and a3=( — 1/2, -\/3/2) represent 
the displacements in a unit cell from the A to B site, from the B to C site, and from the 
C to A site, respectively. In this notation, the first BZ is a hexagon with the corners of 
K=±(2/3)ai, ±(2/3)a2, and ±(2/3)a3. 

In the presence of a magnetic filed (B ^ 0), the hopping term of the tight-binding 
Hamiltonian ([T]) is modified by phase factors from the vector potential A, 

t te'^^i (5) 

where 'jij is the phase factor (the well known Peierls phase) between sites j and i: 

27r f 



= ^ A ■ dl (6) 



with ^Q=h/e being the fiux quantum. 

One can take different gauge potential A mathematically to satisfy the physical confine- 
ment V X A=B. However, the flux through one unit cell is only linearly dependent on the 
external magnetic fleld B, and has no relation with the choice of the gauge potential A. So, 
it is convenient to measure the magnetic fleld in units of the flux quantum per elementary 
triangular plaquette of the kagome lattice. After a simple algebraic calculation, one obtains 
the flux through one triangle is ^=\/3B /A. Thus we can deflne a parameter, called the fllling 
ratio /, as the fraction of a flux quantum throng each triangle, i.e, f=^ /^o=y/3eB / Ah j^. 
The total flux through a unit cell of the kagome lattice is then given as 8/$o- 

When / is a rational number that can be written as f=p/q, where p and q are integers 
with no common factors, the total flux through q triangular plaquette of the kagome lattice 



FIG. 1: (Color online). Schematic picture of the 2D kagome lattice with bond length a. The 
dashed line represents the Wigner-Seitz unit cell, which contains three independent sites {A, B, 
and C). 
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FIG. 2: (Color online). Hofstadter butterfly-like spectrum of the 2D Kagome lattice. The numbers 
recorded in gaps represent the Hall conductance of this system in units of e^/h when the Fermi 
energy lies in these gaps. 

is an integer multiple of $o- In this case, the Hamiltonian ([1]) should be diagonalized in a 
"magnetic" unit cell, which contains q plaquette in order to guarantee k being good quantum 
numbers. As a result, the magnetic Brillouin zone (MBZ) is q times smaller than the usual 
BZ. Furthermore, because of the magnetic translation symmetry, the MBZ has exactly a q- 
fold degeneracy. For convenience we take the gauge potential as A=B(—^{—\/3x + y), |(x — 
■^y),0). Under this gauge, the Peierls phases, which reflect the information about the 
external magnetic field, are only related to n and have nothing with m. The Hamiltonian of 
the kagome lattice now can be diagonalized in the momentum space as the following 3g x 3q 
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matrix, 



9-1 



n=0 



+ hfci^CnC + h^''*ci^iACnC + H.C. 



(7) 



where 



/i„ = t exp(— ik ■ ai) exp 



—iSnf ( n — 
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t exp(— ik ■ exp 



iSnf n + 



12 



and 



/zn*^ = texp(— ik ■ as). 

With the above Hamiltonian ([7]), it is straightforward to obtain the Hofstadter butterfly-hke 
energy spectrum of the kagome lattice, which is shown in Fig. |2j 

From Fig. O one can easily observe three symmetries of this Hofstadter butterfiy-like 
spectrum of the kagome lattice |4i] : (i) The spectrum at / is the same as that at f + j with 
j any integer; (ii) The spectrum is also unchanged on changing / to — /, because if there 
is an eigenstate with energy e for field /, then the corresponding complex conjugate of this 
eigenstate is also an eigenstate with the same energy e for field — /; (iii) The spectrum is 
inverted when the field / changes to / + 1/2. Thus the highest-energy states for field near 
/=l/2 are equivalent to the lowest-energy states near zero field. 

III. THE PROPERTIES OF THE ORBITAL MAGNETIZATION 



To obtain the orbital magnetization of the kagome lattice in an external magnetic field 
and with finite temperature, let us first write down the single-particle free energy as follows 
16|: 



1 



1 + ■ fi„(k) ln[l + e^(^-^"'')l 



(8) 



/5 ~ Jmbz 

Here, fin{k) is the Berry curvature of electronic Bloch states defined by f2.„(k)='i(VkMnk| x 
|VkMnk) with \unk) being the periodic part of Bloch wave for the nth band. Its 



integral over the MBZ gives the topological invariant, namely, the Chern number 
Cn= /mbz d^^QniX). The sum of Chern numbers over the integer occupied bands gives 
the quantized Hall conductance in unit of e^/h (see numbers in Fig. 2). yU in Eq. ([8]) is the 
electron chemical potential, (3=1/ ksT, and e„k is the magnetic band energy. When the field 
varies by an infinitesimal quantity, i.e., from B to B+(5B, then the magnetic band energy can 
be linearly expanded as e„k(B+5B) = e„k(B)— m„k(B) ■ (5B, where m„k is the conventional 
crystal orbital magnetic moment defined by mnk=—i{e/2fi){'VkUnk\ x [H{k) — e„k] l^^Unk) 
in the MBZ. The magnetization is then given by the field derivative at fixed temperature 
and chemical potential, Ai = — {dF/d'B)^rp, with the result 

M = y^ f rf2km„(k)/„(k) (9) 

„ Jmbz 

P „ Jmbz ^ 



Mc + M: 



where /ra(k) is the equilibrium Fermi-Dirac distribution function for nth magnetic Bloch 
band. The first term in Eq. is just a statistical sum of the orbital magnetic moments of 



the carriers originating from the self-rotation of the carrier wave packet 19|, l20l], thus we call 
this term the conventional part of the orbital magnetization. Whereas the second term 
Mn is the Berry-phase correction to the orbital magnetization. This term is of topological 
nature. Interestingly, it is this Berry-phase term that eventually enters the transport current 
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2l|. At zero temperature the general expression ([9]) reduces to [10, 

M = y2 c/'k m„(k) + -a„(k)(/io-e„k) , (10) 
^ 7mbz l " ^ 

where the upper limit means that the integral is over magnetic Bloch states with energies 
below the Fermi energy /io- 

Now with the help of the knowledge on the magnetic energy spectrum obtained from 
diagonalizing the Hamiltonian ([7]), we can numerically obtain the magnetization Ai for 
different field /. First, let us consider the case with field /=l/4. In this case, the energy 
spectrum splits into three bands between which there are two band gaps with the same gap 
width A=1.74 (see Fig. [3](a)). The middle band collapse into a plane in the MBZ. The 
calculated zero-temperature magnetization Ai is plotted in Fig. |3](b). From this figure, one 
can find that with the Fermi energy hq increases from —3.26, the lower band begins to be 




FIG. 3: (Color online), (a) Energy spectrum of the 2D kagome lattice with a field /=l/4. The 
shadow areas are the energy bands. The corresponding magnetization M and its two components 
are respectively drawn in (b) and (c) as functions of the Fermi energy /ig- The shaded areas in 
(b) and (c) are the insulating regions. The dashed and dotted lines in (c) represent and Mq, 
respectively. 

occupied and A4 begins to decrease from 0. When /iq increases to —1.74, the lower band 
is fully occupied and the magnetization decreases to —0.7e/h. When continue increasing 
fit), the system enters the first insulating region. The magnetization then linearly increases 
with the Fermi energy until the middle band is occupied (/io=0~, at this time Ai=1.03e/h). 
When further increasing the Fermi energy, the system immediately comes into the second 
insulating region and the magnetization jumps down to — 1.03e/ h and then linearly increases 
up to 0.7e/h until the upper band begins to be occupied (yUo=l-74). Then the magnetization 
decreases to with the Fermi energy going through the upper band. 



As shown in Refs. 
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18| . the totally different behavior of the magnetization in the 
metallic and insulating regions is due to the different roles and Mq play in these two 
regions. For further illustration, we show in Fig. |3]^c) Mc (dashed line) and Mq (dotted 
line) as functions of the Fermi energy, their sum gives Ai in Fig. |3]^b). One can see that 
overall Mc and Mq have opposite contributions to Ai, which implies that these two parts 
carry opposite-circulating currents. In each insulating regime the conventional term Mc 
keeps a constant, which is due to the fact that the upper limit of the /c-integral of m„(k) is 
invariant as the chemical potential varies in the gap. In the metallic region, however, since 
the occupied states varies with the Fermi energy, thus Mc also varies with /ig, resulting in 
a decreasing slope shown in Fig. [3]^c). The Berry phase term Mq also displays different 
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FIG. 4: (Color online). Energy spectrum and orbital magnetization for /=l/3. 
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FIG. 5: (Color online). Energy spectrum and orbital magnetization for /=l/6. 

behavior between insulating and metallic regions. In the insulating region, Mq linearly 
increases with /iq, as is expected from Eq. (ITOl) . The slope of the Berry-phase correction 
term in insulating regions is proportional to the system's Hall conductance. In the metallic 
region, however, this term sensitively depends on the topological property of the band in 
which the chemical potential is located. On the whole the comparison between Fig. [3]^b) 
and Fig. [3]^c) shows that the metallic behavior of Ai is dominated by its conventional term 
Mc, while in the insulating regime plays a main role in determining the behavior of Ai. 

Now let us consider more complex cases, for example, with a field /=l/3 and with a field 
/=l/6. In these two case, there are more energy bands and gaps appearing in the energy 
spectrum [see Figs, m^a) and[S](a)]. Figures|l](b) andlD^c), and Figs. \5lh) andEI^c), for /=l/3 
and /=l/6, respectively, plot the magnetization Ai and its two components as functions of 
/iQ. Clearly from these figures, one can observe a similar variation of the magnetization by 
changing the electron's fillings (the Fermi energy yUg)- Comparing the case of /=l/3 [Fig. 



S 

6 

-1 



Jo 
|o 

1 


O.DG Q.O^ 






^ f=12/49 

" f=1/4 



-4 -2 2 4 

Fermi energy 



FIG. 6: (Color online). The magnetization Ai as a. function of the Fermi energy hq in the fields 
/=l/4 and /'=12/49. The inset is the enlarged vision of the magnetization A4 at /'=12/49 when 
the Fermi energy varies between 0.05 and 0.1. 
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FIG. 7: (Color online). The two components of the magnetization Ai as functions of the Fermi 
energy //q at /= 12/49. 

m^a)] to the case of /=l/6 [Fig. EJ^a)], one can find that the orbital magnetization in these 
two fields have a relation 

-M(/ = ^,/^o) = -M(/ = ^,-/xo). (11) 

This conclusion is a combined exhibition of the above mentioned symmetries (ii) and (iii) in 
the Hofstadter energy spectrum. 

We have known that the topology of the energy spectrum of the 2D lattice system arises 
from the competition between the periodic potential-induced level broadening and the mag- 
netic field-induced level discretization. Once the external magnetic field changes, the topol- 
ogy of the energy spectrum also changes, which results in the simultaneous prominent vari- 
ation in the magnetization of the 2D system. The above discussions on the magnetization 
for different fields /=l/4, 1/3 and 1/6 are typical examples. Now a question is put forward: 
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If there is a very small perturbation Sf to the field /, then how about the change in the 
magnetization? The answer is the occurrence of fractal structure in the magnetic de Haas- 
van Alphen (dHvA) oscillations, which is very important for the experimental measurement 
of the Fermi surface topology. 

Now we investigate this fractal structure in the magnetic oscillations for the present 
system. For convenience and simplicity, we choose the field /=l/4 and the perturbation 
5/=-l/196 (/'=/ + 5/=12/49). The calculated magnetization M at the fields / and /' are 
plotted in Fig. |6l From Fig. [6] one can observe the following features: (i) In the insulating 
region, the magnitudes of the magnetization Ai have little difference at the fields / and /'. 
The perturbation 6f is more small, more little in the changes of the magnetization Ai; (ii) 
In the metallic regions at the field /, however, there are dHvA oscillations appearing in the 
magnetization Ai at the field / + 6f. The perturbation 6f is more small, the magnetization 
oscillates more rapidly; (iii) The collapsed middle band originally at field /=l/4 now spreads 
at /'= 12/49. And in this middle band there are also dHvA oscillations appearing in the 
orbital magnetization. To see this fractal feature more clearly, we enlarge the magnetization 
in the inset in Fig. |6l (iv) There are similar variations in the two components of At (see Fig. 
[7]). In the insulating regions at /=l/4, there are little difference in the conventional and 
Berry-phase terms. However, both terms exhibit different behaviors in the metallic regions. 
While the Berry-phase term has little difference at different fields as it in the insulating 
regions, there are many quantum steps appearing in the conventional term at /'=12/49. 
The reason for these quantum steps appearing in the conventional magnetization is that the 
topology of the energy spectrum at the field /'=/ + 6f is different from that at the field /. 
When the external field / changes to / + 6f, the eigenbands at the field / split into lots of 
subbands. The perturbation 6f is more small, the number of the splitting subbands is more 
large. Because the upper limit of the fc-integral of m„(k) is invariant as the Fermi level /io 
varies in the subgap, then the conventional term in the subgap keeps a constant, which is 
the reason for the quantum step appearing in the conventional term. On the other hand, 
the integral of the Berry-phase term intimately depends on the Fermi level and has no such 
relation with the energy gaps (or bands). That results in the little difference in the Berry- 
phase term at different fields. So, when one changes the electron's fillings through these 
splitting bands, the total magnetization then exhibits the fractal structure in the magnetic 
dHvA oscillations. 
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FIG. 8: (Color online). The magnetization as a function of the electron chemical potential at 
the field /=12/49 with different temperatures. 

In the above discussions on the orbital magnetization of the 2D kagome lattice, we have 
concentrated on the zero-temperature limit and omitted finite-temperature effect. Now let 
us briefly consider the more realistic cases in which the finite-temperature effect is included. 
Figure [H] plots the orbital magnetization at the field /=12/49 with different temperatures 
kBT=0, 0.005, 0.01, and 0.02, respectively. From Fig. |8l one can find that the magnetic 
oscillations are suppressed by thermal broadening. 

IV. SUMMARY 

In summary, we have theoretically investigated the orbital magnetization of a 2D kagome 
lattice in a perpendicular magnetic field. Here, the orbital magnetization includes a con- 
ventional term and a Berry-phase term, which play different roles in metallic and insulating 
regions. As examples, we have carefully discussed the orbital magnetization and its two 
components at the fields /=l/4, 1/3, and 1/6, respectively. By varying the Fermi energy 
/io, we have obtained the following results: (i) The conventional term and the Berry-phase 
term give the opposite contributions, with their magnitudes increasing (decreasing) with in- 
creasing (decreasing) the Fermi energy in metallic regions; (ii) The conventional term keeps 
unchanged in insulating regions; (iii) The slope of the Berry-phase term in insulating regions 
is proportional to the system's Hall conductance. When the flux is applied near a commen- 
surate one (for example, 1/4), the magnetic dHvA oscillations develop a fractal structure, 
i.e., the orbital magnetization rapidly oscillates when the Fermi energy varies through the 
split subbands. The finite-temperature effect has also been shown to suppress the oscillating 
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amplitude of the orbital magnetization. 
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